Quantum interference and Klein tunneling in graphene heterojunctions 



Andrea F. Young and Philip Kim 

Department of Physics, Columbia University, New York, New York 10027, USA 

(Dated: December 1, 2008) 

The observation of quantum conductance oscillations in mesoscopic systems has traditionally 
required the confinement of the carriers to a phase space of reduced dimensionality^ 1 ^. While 
electron optics such as lensing^ and focusing 6 have been demonstrated experimentally, building a 
collimated electron interferometer in two unconfined dimensions has remained a challenge due to 
the difficulty of creating electrostatic barriers that are sharp on the order of the electron wave- 
length 7 . Here, we report the observation of conductance oscillations in extremely narrow graphene 
heterostructures where a resonant cavity is formed between two electrostatically created bipolar 
junctions. Analysis of the oscillations confirms that p-n junctions have a collimating effect on ballis- 
tically transmitted carriers.—^ The phase shift observed in the conductance fringes at low magnetic 
fields is a signature of the perfect transmission of carriers normally incident on the junctions 10 and 
thus constitutes a direct experimental observation of "Klein Tunneling." 11 



Owing to the suppression of backscatteringi 2 - and its amenability to flexible lithographic manipulation, graphene 
provides an ideal medium to realize the quantum engineering of electron wave functions. The gapless spectrum in 
graphene^ allows the creation of adjacent regions of positive and negative doping, offering an opportunity to study the 
peculiar carrier dynamics of the chiral graphene quasiparticles^iii and a flexible platform for the realization of a vari- 
ety of unconventional electronic devices j 14 i 15 i 16 i 17 Previous experiments on graphene p-n junctio mj 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 
were limited in scope by the diffusive nature of the transport beneath the local electrostatic gates; we overcome such 
limitations by fabricating extremely narrow (~20 nm) local gates strongly capacitively coupled to the graphene channel 
(Fig. la-b). Electrostatics simulations based on finite element analysis (see online supplementary material) show that 
the carrier densities in the locally gated region (LGR) and the 'graphene leads' (GL) — 77-2 and ni, respectively — can 
be controlled independently by applying bias voltages to the top gate (Vtg) and the back gate (Vbg)- The width of 
the LGR, L, is defined as the distance between the two zero density points. As in previous studies 21 , the conductance 
map as a function of Vtg an d Vbg (Fig. lc) can be partitioned into quadrants corresponding to the different signs 
of ni and 77,2, with a lowered conductance observed when 72477,2 < 0. The mean free path in the bulk of the sample, 
l m ^100 nm, was extracted from the relation a = ^-kplm between the conductivity and Fermi momentum, kp. Since 
L < 100 nm within the experimentally accessible density regime, we expect a significant portion of the transport to 
be ballistic in the LGR. 

In the bipolar regime, the diffusive resistance of the LGR is negligible in comparison with the highly resistive 
p-n junctions; as a result, the conductance does not increase with increasing magnitude of the charge density in 
the LGR. 22 We note that the magnitude of this conductance step is only ~ 60% as large as expected for a fully 
ballistic hetero junction even after taking into account the enhancement of the junction transparency due to nonlinear 
screening^ 26 this suggests that there is still a large diffusive component to the transport through the heterojunction. 
Nevertheless, each trace exhibits an oscillating conductance as a function of Vtg when the carriers in the LGR and 
GL have opposite sign. 

The regular structure of these oscillations is apparent when the numerical derivative of the measured conductance is 
plotted as a function of n\ and 77,2 (Fig. 2a). While there is a weak dependence of the oscillation phase on rt\ reflecting 
the influence of the back gate on the heterojunction potential profile, the oscillations are primarily a function of 722, 
confirming their origin in cavity resonances in the LGR. The oscillations, which arise due to interference between 
electron waves in the LGR, are not periodic in any variables due to the strong dependence of the LGR width, L, 
and junction electric field, E, on the device electrostatics. Still, the conductance maxima are separated in density by 
roughly A n2 ^ 1 x 10 12 cm -2 , in agreement with a naive estimate A 77-2 ~ 4 V^ n2 f or the resonant densities in a 
cavity of width L ~ 100 nm. The application of an external magnetic field shifts the phase of the oscillations, with 
individual oscillation extrema moving towards higher density 1 7^2 1 and the transmission resonances appearing to be 
adiabatically connected to the high field Shubnikov-de Haas oscillations (Fig. 2d). 

Graphene heterojunctions offer the opportunity to study an old problem in relativistic quantum mechanics: the 
tunneling of relativistic electrons through a potential barrier^ii In the context of the graphene p-n junction, this 
"Klein tunneling" manifests as the combination of the absence of backscattering with momentum conservation parallel 
to a straight p-n interface: normally incident particles, bound to conserve their transverse momentum, k y = 0, and 
forbidden from scattering directly backwards, are predicted to tunnel through such symmetric potential barriers with 
unit probability. In contrast, particles obliquely incident on a barrier which is smooth on the atomic lattice scale 
encounter classically forbidden regions where the real part of the perpendicular momentum vanishes. These regions, 
which form about the center of individual p-n junctions, transmit obliquely incident carriers only though quantum 
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tunneling, leading to an exponential collimation of ballistic carriers passing through graphene pn junctions, 9 

|T(*g| 2 =e -^v F k 2 y /(eE) ^ 

where vf is the Fermi velocity of graphene. 

Considerable experimental effort has been expended trying to verify eq. (pQ) by matching bulk resistance measure- 
ments across an p-n junctions with their expected value o 19 i 22 i 24 . Such an approach can, at best, provide indirect 
evidence for the theoretically predicted features of chiral tunneling — collimation and perfect transmission at normal 
incidence. In particular, there is no way to distinguish perfect from near perfect transmission from a bulk resistance 
measurement, which is sensitive only to the total transparency of the pn junction. The quantum interference experi- 
ments presented in this letter allow a measurement not only of the magnitude but also the phase of the transmission 
and reflection coefficients. Interestingly, whereas the bulk of conduction in a fully ballistic graphene p-n junction is 
expected to be dominated by normally incident carriers, the absence of backscattering precludes the contribution of 
such trajectories to the Fabry-Perot resonances due to perfect normal transmission at both interfaces. Rather, the 
oscillatory conductance receives its largest contributions from particles incident at angles where neither the transmis- 
sion probability, |T 2 |, nor the reflection probability, \R\ 2 = 1 — |T| 2 , are too large (see, e.g., marker 1 in Fig. 2c). 
Only transmission near such angles contributes to the oscillatory conductance, ensuring the survival of the oscillations 
despite the incident angle averaging and allowing the determination of the width of the angle of acceptance of an 
individual collimating p-n junction. 

In a ballistic hetero junction, the application of a magnetic field bends the carrier trajectories, resulting in an addition 
of an Aharonov-Bohm phase to the interference and a modification of the angle of incidence at each pn junction. As 
was pointed out recently 10 , such cyclotron bending provides a direct experimental signature of reflectionless tunneling, 
which manifests as a phase shift in the transmission resonances of a ballistic, phase coherent, graphene hetero junction 
at finite magnetic field. These resonances are described by the etalon-like ray tracing diagrams shown in Fig 2b. The 
Landauer formula for the oscillating part of the conductance is then 

G osc = e - 2L ' 1 — ^ ^ 2|T+ 1 2 |T_ | 2 1 | cos {Owkb + M rf ) , (2) 

ky 

where T± and R± are the transmission and reflection amplitudes for the classically forbidden regions centered at 
x = ±L/2, Owkb is the semiclassical phase difference accumulated between the junctions by interfering trajectories, 
A9 r f is Klein backreflection phase of the two interfaces, and Ilgr is the mean free path in the locally gated region, 
a fitting parameter which controls the amplitude of the oscillations (see Supplementary Information). 

At zero magnetic field, particles are incident at the same angle on both junctions, and the Landauer sum in eq. 
([2|) is dominated by modes which are neither normal nor highly oblique, as described above. As the magnetic field 
increases, cyclotron bending favors the contribution of modes with k y = 0, which are incident on the junctions at 
angles with the same magnitude but opposite sign (see marker 2 and 3 in Fig. 2c). In the case of perfect transmission 
at zero incident angle, the reflection amplitude changes sign as the sign of the incident angle changes, 10 causing 
a 7r shift in the phase of the reflection amplitudes. Equivalent ly, this effect can be cast in terms of the Berry 
phase: the closed momentum space trajectories of the modes dominating the sum at low field and high k y do not 
enclose the origin, while those at intermediate magnetic fields and k y ~ do (Fig. 2b). Due to the Dirac spectrum 
and its attendant chiral symmetry, there is a topological singularity at the degeneracy point of the band structure, 
k x = k y = 0, which adds a non-trivial Berry phase of ir to trajectories surrounding the origin. As a consequence, the 
quantization condition leading to transmission resonances is different for such trajectories, leading to a phase shift 
in the observed conductance oscillations (i.e., a tt jump in A0 r f) as the phase shifted trajectories begin to dominate 
the Landauer sum in Eq. ([2j). 27,28 For the electrostatics of the devices presented in this letter, the magnetic field at 
which this phase shift is expected to occur is in the range B* ~250-500 mT (see Supplementary Information), in 
agreement with experimental data (see Fig. 3a). As the magnetic field increases further, the ballistic theory predicts 
the disappearance of the Fabry-Perot conductance oscillations as the cyclotron radius, R C1 shrinks below the distance 
between p-n junctions, R c < L, or B~ 2 T for our devices. We attribute the apparent continuation of the oscillations 
to high magnetic field to the onset of disorder mediated Shubnikov-de Haas type oscillations within the LGR. 

In order to analyze the quantum interference contribution to the ballistic transport, we extract the oscillating part 
of the measured conductance by first antisymmetrizing the hetero junction resistance^ with respect to the density 
at the center of the LGR, G~ dd (\ri2\) = G~ 1 (ri2) — G _1 (— 77-2), and then subtracting a background conductance 
obtained by averaging over several oscillation periods in 72,2, G osc = G dd — G Q dd- The resulting fringe pattern shows a 
marked phase shift at low magnetic field in accordance with the presence of the Klein backscattering phase, with two 
different regions — of unshifted and shifted oscillations — separated by the magnetic field B* (see Fig 3a). To perform 
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a quantitative comparison between the measured G osc and eq. ([2]), we first determine the potential profile in the 
hetero junction devices from numerical electrostatics simulations, which information is then input into ([2]) to generate 
a fringe pattern for comparison with experimental data. We choose the free fitting parameter Ilgr = 67 nm for this 
comparison to best fit the oscillation amplitudes. Considering possible degradation of the graphene in and around the 
LGR during the fabrication of the local gates, 29 this value is consistent with the estimate for the bulk mean free path. 
The resulting theoretical calculation exhibits excellent quantitative agreement with the experimental result at both 
zero and finite magnetic field (Fig. 3a-c) both in the magnitude and period of the oscillations. We emphasize that the 
value of L — which largely determines both the phase and amplitude of the oscillations — varies by almost by a factor 
of three over the accessible density range, yet Eq. ([2]) faithfully describes the observed experimental conductance 
modulations in as well as in B. Such remarkable agreement confirms that the observed oscillatory conductance, 
which is controlled both by the applied gate voltage and the magnetic field, results from quantum interference 
phenomena in the graphene heterojunction. Moreover, the oscillations exhibit a phase shift at B* ~ 0.3 T which is 
the hallmark of perfect transmission at normal incidence, thus providing direct experimental evidence for the "Klein 
tunneling" of relativistic fermions through a potential barrier. 

Finally, we turn our attention to the temperature dependence of the quantum coherence effects described in the 
text, which we observe at temperatures as high as 60 K (Fig. 3d). An elementary energy scale analysis suggests that 
the phase coherence phenomena should be visible at temperatures of order ~ 100 K, when thermal fluctuations 
become comparable to the phase difference between interfering paths. In addition, the oscillation amplitude is sensitive 
to the carrier mean free path, and we attribute the steady waning of the oscillations with temperature to a combination 
of thermal fluctuations and further diminution of the mean free path by thermally activated scattering. The mean 
free path in clean graphene samples can be as large as ~1 /im^ and a reduction of the width of the heterostructure 
L by an order of magnitude is well within the reach of modern fabrication techniques; consequently, technological 
improvements in the fabrication of graphene hetero junctions should lead to the observation and control of quantum 
coherent phenomena at much high temperatures, a crucial requirement for realistic, room temperature quantum device 
applications. 

Methods 

Graphene sheets were prepared by mechanical exfoliation 13 on Si wafers covered in 290 nm thermally grown Si02- 
Ti/Au contacts 5 nm/35 nm thick were deposited using standard electron beam lithography, and local gates sub- 
sequently applied using a thin (~ 10 nm) layer of hydrogen silsesquioxane (HSQ) as an adhesion layer— for low- 
temperature atomic layer deposition of 20 nm of Hf02, a high-A: dielectric (e ~ 12) (see Fig. lb). Palladium top 
gates not exceeding 20 nm in width were deposited in order to ensure that a sizeable fraction of conduction electrons 
remained ballistic through the LGR. Leakage current was measured to be < 100 pA up to Vtg = ±15 V. All data 
except Fig. 3d was taken from the device depicted in Fig la, which had a measured mobility ~ 5,000 cm 2 /V sec. 
Fig. 3d was taken from a similar device in a four terminal Hall bar geometry; additional data from this device is shown 
in the supplementary materials. Several other similar devices were also measured, showing qualitatively similar be- 
havior. The conductance of the graphene devices was measured in a liquid helium flow cryostat at 4.2- 100 K using a 
standard lock-in technique with a current bias of .1-1 fiA rms at 17.7 Hz. Unless otherwise specified, all measurements 
were done at 4.2 K. The ratio Ctg/Cbg ~ 12.8 was determined from the slope of the Dirac ridge with respect 
to the applied voltages, and similar values were obtained from the analysis of the period of the Shubnikov-de Haas 
oscillations in magnetic field, which also served to confirm the single layer character of the devices. Finite element 
electrostatics simulations were carried out for the measured device geometries described above with the thickness and 
dielectric constant of the HSQ adjusted such that the simulations matched the observed values of Ctg/Cbg- The 
shape of the potential and the strength of the electric field E used in fitting the experimental data were constrained 
to lie within the confidence interval of the simulations, which in turn were largely determined by uncertainty in the 
device geometry. 
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FIG. 1: Graphene heteroj unction device schematics and conductance measurements, a False color scanning electron microscope 
image of a typical graphene heteroj unction device. Electrodes, graphene, and top gates are represented by yellow, purple and 
cyan, respectively. The scale bar is 2 /im. Inset: high magnification view of top gate. The scale bar is 20 nm. b Schematic 
diagram of the device geometry. The electrostatic potential created by the applied gate voltages, Vbg and Vtg, can create a 
graphene heteroj unction of width L bounded by two p-n junctions, c The inset shows the conductance as a function of Vtg 
and Vbg- The main panels show cuts through this color map in the regions indicated by the dotted lines in the inset, showing 
the conductance as a function of Vtg at fixed Vbg- Traces are separated by step in Vbg of 1 V, starting from ±80 with traces 
taken at integer multiples of 5 V in black for emphasis. 



6 




FIG. 2: a dG/dri2 as a function of m and 712. Scale bar is in arbitrary units, b Schematic of trajectories contributing to quantum 
oscillations in real and momentum space. The dominant modes at low magnetic field (marker 1) give way, with increasing B, 
to phase shifted modes with negative reflection amplitude due to the inclusion of the non-trivial Berry phase (marker 3), near 
k y = 0. The original finite k y modes are not yet phase shifted at this field (marker 2), but due to collimation, they no longer 
contribute to the oscillatory conductance, c The prefactor in the Landauer sum, |T+| 2 |T_ \ 2 R+R- — |T+| 2 |T_| 2 |i?+||i?-|e 2A6,r/ , 
as a function of B and k y , plotted for the experimental parameters at Vbg=50V for 712 — 3.5 x 10 12 cm -2 . The markers 
correspond to the trajectories shown in b. Regions of negative sign correspond to trajectories containing the Klein backscattering 
phase shift, d Main panel: Magnetic field dependence of dG/dri2 at Vbg = 50 V. Inset: Similar data taken at Vbg = —50 V. 
The magnetic phase is proportional to the sign of the carriers; as a result, the oscillation extrema precess in opposite directions 
for opposite signs of carriers in the LGR. 
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FIG. 3: Comparison of experimental data with the theoreticail model, and temperature dependence, a Magnetic field and 
density dependence of the oscillating part of the conductance at Vbg = 50 V. G osc as extracted from the experimental data 
(left panel) shows good agreement with a theoretical model accounting for nonlinear screening 26 (right panel) over a wide range 
of densities and magnetic fields, b Cuts taken at B = 0, 200, 400, 600, and 800 mT, corresponding to the colored dashed 
lines in a; the dots represent data, the smooth lines the result of the simulations. The sudden phase shift that signals the 
presence of perfect transmission is indicated by dashed arrows. Curves are offset for clarity, c Temperature dependence of the 
oscillation amplitude in a similar device to that presented in the rest of the text. Main panel: The oscillations (different curves 
are offset for clarity) weaken with rising temperature, and are not observed above 80 K. At 4 K, the conductance modulations 
contain both the ballistic oscillations as well as aperiodic modulations due to universal conductance fluctuations which quickly 
disappear with increasing temperature. Inset: Averaged amplitude of several oscillations, normalized by the amplitude at 
T = 4 K. 
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Supplementary Information 
Details of theoretical model 

The Landuaer formula for the oscillating part of the conductance is obtained from the ray tracing scheme shown 
in Figure 2a of the main text. The transmission amplitude through the entire junction is taken to be a product 
of the transmission amplitudes at the two interfaces with a phase factor corresponding to the semiclassical phase 
accumulated between the junctions. This semiclassical phase difference for neighboring trajectories is 



J^J^\n(x)\-(k y -^Bx)\ (3) 



OwKB = 3? 

-L/2 

where we take the real part to account for the fact that in general the classical turning points shift from their values 
of ±L/2 defined for k y = 0, B = 0. In addition, there is a nonanalytic part of the phase associated with the the 
vanishing of the reflection coefficients at opposite interfaces^ which can be nontrivial at finite magnetic field, 

where H(x) is the step function (H(x) = 1 for x > 0, H(x) = otherwise). This phase jump is equivalent to a sign 
change in the reflection coefficient as the incidence angle crosses zero, and implies that the transmission probability 
at normal incidence is unity. 

At a p-n junction, assumed to be smooth on the lattice scale, the transmission amplitude is exponentially peaked 
about normal incidence. 9 The principal effect of a weak magnetic field on the transmission through a single p-n 
junction is to modify the incident angle at the two junctions due to cyclotron bending of the trajectories. Choosing 
the Landau gauge where A = Bxy with k y the conserved transverse momentum in the center of the junction (x = 0), 
the transmission and reflection amplitudes at the junctions located at x = ±L/2 are^ 

/ \ 2\ 

ttHvf ( 7 eBL\ 



T± = exp — -±- [ky ± — - R± = exp (iirH (-ky T eBL/(2h))) ~ \T±\ 2 (5) 



^ 2eE \ Ky 2h ) / 

Defining the total phase 6 = Owkb + A#, and taking into account the damping, due to scattering, of the particle 
propagators between the junctions, we can write the Landauer conductance of the hetero junction via the canonical 
Fabry- Perot formula as 



Ap 2 „ T 1 , t P -L/(2i LGR ) 

4t x i+i-e (6) 



1 - \R + \\R-\e i9 e- L / l ^ R 



It is difficult to separate diffusive from ballistic effects in the bulk conductance. However, the contribution of 
diffusive effects to quantum interference effects are strongly suppressed at B = 0, rendering useful the definition 
Gqsc = G — G, where G denotes the conductance averaged over the accumulated phase. Multiple reflections are 
suppressed both by the finite mean free path as well as the collimating nature of the junctions; as |T±|, \R±\ < 1, 
higher order products of transmission and reflection coefficients are necessarily decreasing. Utilizing this fact, we 
can expand the denominator in eq. ([6]) and subtract all nonoscillating terms to get the leading contribution to the 
quantum interference, 



°osc = ^-^2\T+\ 2 \T-\ 2 \R+\\R-\cos(6) e - 2L ' 1 —. (7) 

k 

K y 

This formula does not take into account inhomogeneities in the applied local gate potential due to the uneven 
width of the top gate or thickness and crystallinity of the dielectric. Such disorder can lead to significant damping of 
the oscillations even in a completely ballistic sample^, and makes the estimate for the mean free path derived from 
matching the observed amplitude of the oscillations a lower bound. Nevertheless, that the mean free path in the LGR 
is shorter than that in the GL is consistent with recent experiments (unpublished) which find a strong enhancement 
of the Raman D-band in graphene after it is covered in HSQ and irradiated with electrons at energies and doses 
comparable to those used during top gate fabrication^ 
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Electrostatics simulations 

The theoretical model described in the previous section takes as an input the potential profile across the hetero- 
j unction. In order to perform a quantitative comparison between experiment and theory, we first determine these 
parameters by numerical simulation of our device electrostatics using Comsol Multiphysics, a commercial finite el- 
ement simulation software package. Scanning electron microscopy images of the device chosen for the quantitative 
comparison show the top gates to be ~ 20 nm wide, while the capacitive coupling of the top gate is found to be 
Ctg ~ 12.8Cbg- Graphene is treated as a perfect conductor covered by 10-20 nm of hydrogen silsesquioxane (HSQ), 
with dielectric constant e ~ 2 — 5, and 20 nm of Hf02, e ~ 10 — 15. We choose the dielectric constant of the HSQ 
in each simulation to ensure a match between the calculated and measured capacitive coupling, determined as the 
relation between the applied top gate voltage and the maximal density reached in the LGR. The dielectric constant 
of atomic layer deposition (ALD) grown Hf02 layer was consistent with a separate capacitance measurement on a 
similarly prepared thin film. For this range of sample parameters, we find the density profile in the hetero junction to 
be well approximated by 

n(x)= CT ^ 25 +C BG V BG , (8) 
1 + \x/w\ 

where w ~ 45-47 nm is the effective width of the potential and the gate potentials Vtg and Vbg are coupled to the 
charge density through the capacitances Ctg ~ 1490 aF/im 2- and Cbg ~116 aF/im -2 . Here the exponent of x/w 
is chosen empirically, and produces a potential profile that deviates by less than 10% from that produced by the 
electrostatics simulations. 

Fig. 4a shows density profiles in the graphene hetero junction which produce, via the Landuaer formula described 
in the previous section, conductance fringes matching those observed in experiments done at Vbg — 50 V. From these 
density profiles, the potential profile, distance between classical turning points, and density gradient at the center of 
the p-n junctions is determined. We note that the distance between pn junctions, which determines the boundary of 
integration in equation ([3|), changes by a factor of two over the experimentally accessed density range, making the 
excellent quantitative fits obtained for the oscillation phase particularly convincing. 
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FIG. 4: Density profiles for Vbg = 50 V over a range of top gate voltages (from eq. [5}. Top gate voltage in 1 V increments. 
The width of the central region (L) depends strongly on applied voltage. 
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Quantitative estimate of the electric field at the p-n interfaces 

In this section, we compare the measured data (Fig. 5a) with numerical simulations to extract the strength of the 
electric field in an individual pn junction. It follows from robust semiclassical arguments 9 that, in the experimentally 
realized situation of p-n junctions smooth on the scale of the lattice constant, the collimation at an individual junction 
should be a Gaussian function of k y (see eq Ej). The important parameter is the electric field, E, in the junctions, 
which is given, after taking into account the absence of linear screening in graphene near the charge neutrality point, 
by 26 

eE = 2.1hv F n ,2/s . (9) 

where n' is the density gradient across the junction. As is evident from the simulations (Fig. 5b), the nonlinear 
screening correction to the electric field gives a better fit to the experimental data than either the non-exponential 
collimation produced by atomically sharp barriers — which appears to contain higher order resonances — or the weaker 
field that results from neglecting nonlinear screening near the Dirac point. 

To make this comparison more quantitative, we perform several simulations in which the nonlinear screening result, 
([9]), is scaled by some prefactor, 77. For a rough comparison, the linear screening corresponds to 77 < .5, while the 
step potential corresponds to 77 ^> 20. As explained in the main text, the magnetic field dependence contains an 
abrupt phase shift at finite magnetic field as the finite k y modes cease to contribute to the oscillations and k y = 
modes — which carry with them the additional Berry phase of tt — become the dominant contribution to the oscillatory 
conductance. With increasing magnetic field, a fully ballistic model predicts the gradual ebbing of these phase shifted 
oscillations as the cyclotron radius becomes comparable to the junction size. The field at which the 7r-phase shift 
manifests is tied to the degree of collimation of the transmission at each p-n interface. Because this phase shift is 
rather abrupt, we can define the transition magnetic field, as the field at which the values of the the oscillation 

prefactor |T+| 2 |T_| 2 |i?+||i?_| for zero and finite k y modes become comparable, giving B* oc yj ' e 2^ L 2 - Since B* 

depends strongly on the junction electric field, it allows us to extract this field from the experimentally observed 
oscillation phase shift. Defining B* as the field at which G OS c(n2,B = B*) = for fixed density ri2 such that 
G osc (ri2,B = 0) is an extremum, we can estimate 77 = .9 ± .3. 

In accordance with the ballistic theory, the oscillations peak at zero magnetic field, and then have a second maximum 
after the phase shift at finite magnetic field. The relative height of these two maxima can be used to estimate the 
electric field E. Higher collimation suppresses the contribution of the modes near k y = at finite £?, since this 
feature is entirely generated by modes not normally incident at either interface. Higher collimation thus corresponds 
to an effectively more one-dimensional channel for interference effects, leading to the more effective destruction of the 
oscillations by the Lorentz force, which serves to push the particles out of the narrower acceptance angles at each 
junction. By taking the average value of Gg ax jG^ ax and comparing it with the simulations for a variety of values of 
77, we can extract 77 «1.6±.3. We note that the apparent discrepancy between this and the method previously outlined 
is likely due to the contribution of momentum non-conserving processes to the oscillations at finite leading to an 
overestimate of 77. While neither of the methods above is free of systematic errors, they confirm the importance of 
nonlinear screening in determining transport through graphene p-n junctions. 
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FIG. 5: a Top panel: oscillatory conductance as a function of ri2 and magnetic field at Vbg=50V. b The oscillation prefactor 
\T + \ 2 \T-\ 2 R+R- (column 1) and resulting oscillatory conductance (column 2) as a function of magnetic field for a variety of 
collimation models. The best fit to the data is achieved by accounting for the nonlinear screening (top panel); the simulations 
resulting from naively linearizing the potential between the extrema (middle panel) and considering the algebraic collimation 
resulting from a step potential (lower panel; note the different scale in left) show features incompatible with the observed data. 
The width of the central region is adjusted to be 46-48 nm in the simulations in order to match the phase of the zero field 
oscillations. The amplitude is fit by setting the mean free path in the Landauer formula to be 67 nm in the top panel, 60 nm 
in the central panel, and 300 nm for in the lower panel, c Transmission probability as a function of angle at zero magnetic 
field for the pn junctions with (red, solid) and without (green, dotted) nonlinear screening, and for the step potential (black, 
dashed) . 



Disordered heterojunctions 

Several measured devices showed Fabry-Perot resonances; data from a second device (the same used in Fig. 3d of 
the main text), is shown in Fig 6. Most of features discussed in text are present, including an oscillatory conductance 
that can be tuned be magnetic field. In this device, however, as in the majority of devices, the Fabry-Perot resonances 
appear to be mixed with other, irregular, conductance fluctuations. This behavior is most evident in the magnetic 
field dependence of the oscillations (Fig. 6c): although the phase shift is still evident in some of the transmission 
resonances, several oscillation periods appear to be intermixed, and there is no adiabatic crossover to the Shubnikov-de 
Haas oscillations as observed in the high quality device discussed in the main text. In addition, we discovered that 
the more disordered heterojunctions devices exhibit transmission resonances even when the overall doping of the LGR 
was of the same sign as that in the GL (Fig. 6b). We interpret these effects as a combination of higher disorder 
concentration between the p-n junctions causing enhanced universal conductance fluctuations and inhomogeneous 
gate coupling, which has the effect of causing an averaging over several Fabry Perot fringes. 10 
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FIG. 6: Data from a second device, a Inset: Conductance, measured in the four wire geometry, as a function of top and back 
gate voltages. Main Panel: Derivative of conductance with respect to top gate voltage (arbitrary units), b Oscillations in the 
unipolar regime. The oscillations are visible in the derivative of the conductance with respect to top gate voltage, and as kinks 
in the steeply falling background conductance as the central region approaches charge neutrality, c The irregular character of 
the oscillations becomes apparent in the magnetic field dependence of the oscillation extrema. Although the phase shift is still 
visible, several periods seem intermixed, and there is no apparent adiabatic connection to the Shubnikov-de Haas oscillations. 



